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Abstract 



In this paper we propose a novel algebraic and geometric description 
for the dissipative dynamics. Our formulation bears some similarity with 
the Poisson structure for non-dissipative systems. We develop a canonical 
description for constrained dissipative systems through an extension of the 
Dirac brackets concept, and we present a new formula for calculating Dirac 
brackets. This formula is particularly useful in description of dynamical 
systems with many second-class constraints. After presenting necessary for- 
mal background we illustrate our method on several examples taken from 
particle dynamics, continuum media physics and wave mechanics. 



1 Introduction 



The first systematic attempt to provide mathematically consistent quanti- 
zation procedure for constrained systems was given by P.M. Dirac who 
derived a formal "replacement" for the canonical Poisson brackets which 
plays today a fundamental role in the canonical formalism for constrained 
Hamiltonian systems on both classical and quantum levels. In spite of the 
considerable attention paid to that formula in the mathematical literature 
1^, ^, 0] , and several attempts to use the Dirac brackets in the quantization 
of gauge invariant systems, eg. [||, ^ etc., until recently there were few at- 
tempts to actually use this formalizm in more conventional applications. We 
have recently provided few examples of these applications in classical and 
continuum mechanics |^. 

The canonical formalizm applies to the conservative systems. These sys- 
tems form a relatively small sub-class of interesting physical systems since 
most of the others describing phenomena at some effective rather than fun- 
damental level are dissipative. In the past several attempts were done to 
describe dissipative classical mechanics in a fashion similar to its canonical 
description. One of these attempts, so called metriplectic approach |^, ^] 
was advocated as a natural extension of the mixed canonical- dissipative dy- 
namic proposed by Enz In this paper we extend further the concept of 
metriplectic dynamics to what we shall call semimetric-Poissonian dynam- 
ics, which is nothing but a natural combination of semimetric dynamics (a 
dissipative part) and Poissonian dynamics (a conservative part). We pro- 
pose a canonical description for constrained dissipative systems through an 
extension of the concept of Dirac brackets |||] developed originally for conser- 
vative constrained Hamiltonian dynamics, to the non-Hamiltonian, namely 
metric and mixed metriplectic, constrained dynamics. It turns out that this 
generalized unified formula for the Dirac brackets is very useful in the de- 
scription and analysis of a wider class of dynamical systems. To proceed 
with our approach we develop a new formula for calculating Dirac brackets 
which is particularly effective in finding equations of motion and constants 
of motions for systems with many constraints. 

In order to make this paper self-contained we include in Section II a 
short "primer " to the Poisson geometry, dynamics and the Dirac brackets. 

The rest of the paper is organized into four sections. In Section III we 
discuss semimetric dynamics and some elementary features of the related 
mathematical structures like semimetric algebras and SJ-identity the late 
should be regarded as dynamical symmetric version of the Jacobi identity. 
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We develop symmetric concepts in the analogy with those in the Poisson 
category in Appendix. 

In Section IV we discuss semimetric dynamics subject to some con- 
straints. We derive a symmetric analogue of the Dirac brackets and pro- 
vide its geometric interpretation as induced metric on some submanifold of 
a Riemannian manifold. We also present a new effective algorithm for cal- 
culation of the Dirac brackets in both symmetric and antisymmetric cases. 
Few examples, for finite and infinite dimensional cases, are also discussed in 
some detail. 

In the section V we discuss semimetric-Poissonian dynamics which is a 
combination of semimetric dynamics, discussed in the section III, and Pois- 
sonian dynamics mentioned in the section II. This section also includes the 
extended Dirac approach to constrained semimetric-Poissonian which is a 
combination of the Dirac approach to constrained semimetric dynamics dis- 
cussed in the section IV and the usual to constrained Poissonian dynamics. 

In Section VI we discuss interesting physical examples: the dissipative 
formulation of the Schrodinger equation due to Gisin the Landau- 

Lifshitz-Gilbert equations for damped spins ||8|, dynamics of damped 
rigid body - including among other the Morrison equation of damped rigid 
body We also discuss a novel description of the incompressible viscous 
fluid. 

Some important but not crucial mathematical aspects of the extended 
canonical formalism for constrained metric and mixed semimetric-Poissonian 
systems will be given in the forthcoming publication [13|. Computational 
aspects of the Dirac brackets - symmetric, antisymmetric or mixed - will 
be publish shortly [0]. 



2 Poisson geometry, dynamics and Dirac brackets 

From the algebraic point of view the Poisson algebra is a linear space 
equipped with two structures: 

i) the commutative algebra structure with the (associate) multiplication 

T ^ T ^ T ^ the product of two elements /, g is denoted simply by fg^ 

ii) the Lie algebra structure with the Lie bracket {•,•}: .7^ x ^ ^ .F, the 
product of two elements /, g is denoted by {/, 5}, 

related to each other by the Leibniz rule: 

{/9,M = /{5,M + {/,M5- (2.1) 
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The Lie bracket of a Poisson algebra is called the Poisson bracket. 

The Poisson manifold is a smooth manifold M for which the commuta- 
tive algebra of smooth functions on M, C°^{M), is equipped with the Pois- 
son bracket. The Poisson bracket {•, •} acts on each function as a derivation, 
thus there exists a contravariant (2, 0)-tensor 11 such that {/, g} = Il{df, dg) 
for every functions f,g. In the local coordinates {z^) 

m r. 

U,g}{z) = Y.^'\z) dj d,g ,dk^-^. (2.2) 

The tensor 11 which defines a Poisson bracket is called a Poisson tensor. 
The antisymmetry of Poisson bracket implies that tensor 11 must be anti- 
symmetric, so n*^ = — n-'*. The Jacobi identity requires that 

N 

* diW'' + ^' diU^' + U^^ diW^ = . (2.3) 

1=1 

A Hamiltonian vector field generated by a function /t is a vector field defined 
by Xh{f) = {/, h} for V/. All flows generated by Hamiltonian vector fields- 
Hamiltonian flows, preserve the Poisson structure. 

Poisson dynamics or generalized Hamiltonian dynamics is a dynamics 
generated by some Hamiltonian vector field, for which the Hamiltonian func- 
tion plays physically specific role. 

Let TM and T*M denotes the tangent, and cotangent bundle of the 
manifold M. The Poisson tensor H induces a bundle map ntt : T*M TM 
which is defined by Ii^{df ) := Xf for all functions /. The rank of a Poisson 
structure at point z is defined to be the rank of n| : T*M — > T^M, which is 
equal to the rank of matrix IT^^z) in the local coordinates (z^)- A Poisson 
structure with constant rank equal to the dimension of the manifold M is 
called symplectic or nondegenerate. In this case the "inverse map" of H, 
denoted by u, is a symplectic 2-form and uj{Xf,Xg) = {f,g} = Il{df,dg). 
The Darboux theorem states that for every symplectic structure there ex- 
ists, locally, a canonical coordinates system (xi, . . . , Xk,Pi, . . . ,Pk) such that 
H = ^ A ^, (iu = dx A dp), or equivalently, {xi,Xj} = {pi,Pj} = 0, 

The invariance of the Poisson structure under Hamiltonian flows implies the 
constancy the tensor H rank along the orbits of such flows. The orbit of 
each point of M under the action of all Hamiltonian flows forms a symplec- 
tic manifold called a symplectic leave. Since M is a union of such orbits 
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every Poisson manifold is a smooth union of disjoint connected symplectic 
manifold (symplectic leaves) of various ranks. 

The splitting theorem [jl5| for Poisson manifolds, states that locally every 
Poisson manifold is the product of a symplectic manifold and a Poisson man- 
ifold with zero rank. In the other words, locally in the neighborhood of the 
point zq there always exist canonical coordinates system: 

(xi, . . . ,Xk,Pi, . . . ,Pk,Z2k+i, ■■■,Zn) such that {xi,Xj} = {pi,Pj} = 0, {xi,pj} = 
5ij, {xi, zi} = {pj,zi} = 0, {zr, Zs] = Ars and As(zo) = 0. 
The map : Mi M2 between two Poisson manifolds is called Poisson 
mapping iff {f o (p,g o ip}i = {f,g}2 o 99. The Poisson mapping is a natu- 
ral generalization of the well-known from classical mechanics notion of the 
canonical transformation. 

In the usual formulation of the classical mechanics the constrained dynamics 
can be visualized geometrically as the dynamics on some submanifold of the 
system phase space. Similarly, the constrained Poisson dynamics can be rep- 
resented as such on some submanifold of the Poisson manifold. However, it 
is not always possible to define induced Poisson structure on a submanifold 
and therefore we have no obvious way how to generalize constrained Pois- 
son dynamics. If the Poisson structure is non-degenerate (symplectic case) 
then on each submanifold there exists an induced two form which becomes 
symplectic if it is non-degenerate. Further for an arbitrary submanifold N 
of a symplectic manifold (M, co) there always exists a maximal submanifold 
N' C N such uj^^i is non-degenerate, then A^' has an induced symplectic 
hence Poisson structure. 

Dirac IQ has proposed an algebraic procedure to deal with constrained dy- 
namics. First consider any finite set of (linearly independent) constraints 
A = {ai, . . . ,ak} <Z T = C°°{M) and define weakly- vanishing (or weakly- 
zero) elements as linear combinations of constraints with arbitrary functions 
coefficients, i.e. I = {X^i/j^i; where fi G J^,ai € A}. 
Element f £ J- is called first-class with respect to the set of constraints A 
iff it has weekly-zero bracket with all constraints, i.e. \/a A : {/, a} € X. 
Otherwise, it is called second-class. The set of all first-class elements de- 
noted by J^i {A) , forms a linear subspace of J- and the set of all second-class 
elements, denoted by J^2iA). 

This classification divides the set of constraints A into two subsets: first-class 
constraints Ai = Ad Ti (A) and second-class constraints A2 = Ad T2 (A) . 
The number of second-class constraints must be even A2 = {©i, . . . , 02s}- 
Dirac has proven the Gramm matrix of second-class constraints, [{Qi, Qj}] = 
W, is weakly non-degenerate. This allowed him to define new antisymmetric 
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Leibniz bracket, known as the Dirac bracket: 

2s 

{f,9}D = {f,g}-J2if'®^}Ciji®j^9} (2.4) 

where C = [Cij] = W"^ is an inverse matrix of W. Using ( p.4D one can 
check that the Dirac brackets posses ah the required properties of the Poisson 
brackets. The (algebraic) proof of the Jacobi identity is difficult. One can 
easily check that all the second-class constraints are Casimirs with respect 
to the Dirac bracket, i.e. {©fc,/}D = for all /. 

3 Semimetric manifolds and semimetric dynami- 
cal systems 

In this section we introduce a concept of semimetric algebras and semimetric 
manifolds which play a similar role in the description of dissipative systems 
as the Poisson algebras and the Poisson manifolds in the description of con- 
servative dynamics. In the last two sections of this work we will show that 
the semimetric structure together with the Poisson structure are sufficient 
for "canonical" description of a wide class of dissipative dynamical systems. 

Let X be non-empty set. A semimetric bracket on the linear space of 
real functions defined on X, namely = Fun(X), is a bilinear operation 
-< •, • T y. T ^ T which satisfies the following requirements: 

i) It is symmetric: V/, (/ € : ^ /, g = < g, f 

ii) It satisfies the Leibniz rule: \/f,g,h € : -< fg,h >-=^ f,h >- 
g + f^g,hy. 

ill) It is non-negative definite: V/ ^ T : -< f,f ^> 0, i.e. the function 
-< f, f y is non-negative definite function, Vx G X -.^ f, f y (x) > 0. 

Note that if a bilinear operation on satisfies conditions i), ii) it is called 
pseudo-metric bracket or symmetric Leibniz bracket. A symmetric Leibniz 
bracket which satisfies the condition 

iii*) Positive definite: £ T : ^ /, / ^> and ^ /, / ^= iff 
/ = Const (at least locally), 
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is called metric bracket. 

Symmetric Leibniz algebra is a linear space J-" equipped with two struc- 
tures: commutative algebra structure with the (associate, commutative) 
multiplication and symmetric (non-associate, non-commutative) structure 
and these two structures are related by the Leibniz rule. Semimetric algebra 
is a symmetric Leibniz algebra whose bracket is semimetric. 
As we shall see the semimetric algebra can be used to describe the (wide 
class) of the dissipative classical systems akin to the description of the non- 
dissipative dynamics by means of the Poisson algebra. 

From now we assume that X is a smooth finite dimensional manifold, namely 
Phase space, and = C°°{X) is a space of all smooth functions on X. 
There is a correspondency one to one between the symmetric Leibniz brack- 
ets on the space of functions and the symmetric tensors on the manifold X: 
-< f,g y = G{df, dg) where G is a contravariant tensor field of the type (2, 0) 
on X. In the local coordinates {z^), each symmetric tensor is of the form 
G(z) = X^ij G'*"' (^) ^ 'SD ^ , where G*-' = G-'*, hence each symmetric Leibniz 
bracket locally must be of the form: 

^/,9^(^) = jl G^^{z)^^^ , f,geC^{X). (3.1) 
«,i=i 

A symmetric bracket becomes semimetric bracket iff it is non-negative, i.e. 
the matrix [G*-'] is non-negative definite. 

A semimetric bracket is called a metric bracket if G is positive definite (non- 
negative and non-degenerate), i.e. with constant maximal rank, rank G = 
dim X. In the general the tensor G may have a nonconstant rank which 
depends on points. 

Definition 3.1 Semimetric manifold is a smooth manifold M for which 
the commutative algebra of smooth functions on M is a semimetric algebra, 
i.e. it is equipped with a semimetric bracket. Geometrically, semimetric 
manifold can be viewed as a pair (M , G) where G is semimetric tensor (or 
cometric tensor), i.e. symmetric, non-negative definite: G{df,df) > for 
every function /, contravariant (2,0)-tensor field. 

Example 3.1 Tensor G(z) = J2i=i £-^£--^]=k+i A]® A- symmetric 
and non-degenerate, but it is not non-negative definite. Tensor G(z) = 

^Note that it is true for smooth functions C°°{X), but false for C'°(X). 
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Eil J- ® J- for fc < n, and tensor G(z) =zf^<S>^ + (zizi)^ ® + 
J27=3 '5^i ^'^^ non-negative definite, but degenerate. □ 

A dissipative vector field generated by the function his a vector field defined 
by Xj^{f) f, h >~ for all functions /. The flow generated by dissipative 
vector field should be called a dissipative flow. Dissipative flows essentially 
differ from the Hamiltonian counterpart: they do not preserve the symmet- 
ric structure. 

If the (2, 0)-tensor G is positive-definite, there exists symmetric (0, 2)-tensor 
Q, its inverse, such that G{df, dh) = Q{Xj^ , Xj^), which is exactly a Rieman- 

nian metric tensor. In the local coordinates (zk), if: G{z) = J2ij ^*''(^)^® 

the tensor G{z) = Y,ij Gij{z)dzi ® dzj where J2j G'^^Gjk = SI- 
The concept of semimetric manifold then is a natural generalization concept 
of Riemann manifold. It is analogous to a generalization from symplectic 
manifold to the Poisson manifold. 

Similarly as in the Poissonian case, it is not always possible to define in- 
duced semimetric structure on a submanifold, hence there is no obvious way 
to describe canonically constrained semimetric dynamics. 

The map F : Mi M2 between two semimetric manifolds is called 
semimetric mapping iff it maps the semimetric structures, i.e. -< f o F,g o 
F yi=^ f,g ^2 oF. 

Proposition 3.1 Let {!F, •, -< •, • y) be semimetric algebra, 
a) The Schwartz inequality holds 

yf,geJ^: -<f,fy^g,gy><f,gy^. (3.2) 



b) Let /i, /2, • • • , /n be arbitrary elements of J^. Then the square matrix 



Gram{fi, . . . , fn) = 



-< fiJi y 
-< /2,/i y 

~< fn, fl >- 



-< fiJn y 

-< /2,/n >- 
~^ fm fn y 



(3.3) 



is non- negative definite. In particular, det Gram{fi, . . . , fn) > 0. Fur- 
thermore, if ^ is a metric algebra, then det Gram{fi, . . . , fn) = iff 
{fi}i=i are affine linear dependent. 
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Proof. a) Indeed, for each real number A G R the expression <^ 
/ - \gj - Xg y= A2 ^ g,g > 2A ^ f,g ^ + ^ fj >- is non- 
negative quadratic form in the real number A. Hence the discriminant 
A = 4 ^ f,g y'^ —4 -< /, / g,g y< must be non-negative. 
If T is metric algebra, then -< f,g ^^=-< /, / g,g y iE f,g are afhne 
linear dependent, i.e. f — Xg = Const. 

b) For each vector a = (oi, 02, . . . , a„) G R", denote a • f = CLifi, we 
have a'^ Gram{fi, . . . , /„) a =-< a • f , a • f ^> 0. 

Q.E.D. 



Example 3.2 

1. The natural Euclidean metric of n-dim Euclidean space X = R" in- 
duces a natural metric structure on C°°{X) 

^/,<7^(x) = X]^|^, where /, g G C-(R"), x G R" . (3.4) 

i=l * * 

We shall call the metric bracket ( p. 4]) the Euclidean metric bracket. 

2. Here is a simple but general construction of the semimetric structure 
in the space of smooth functionals over a Hilbert space. Let X be 
a Hilbert space with a scalar product < • | • > and A be a linear 
operator on X. One can define a semimetric structure on the space of 
all smooth functionals over X as following 

For instance, let X be a Hilbert space of functions X = {i : R" — > R*^} 
with a scalar product 




f = (/i,/2,---,/d), g = {gi,92,---:9d) ex. (3.6) 
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The above scalar product on X defines a semnnetric bracket on the 
space of smooth functional over X by 



~-^^J dfi{x) 6fj{x) 

where G C°°(X). 



To be specific let d = n and G be a differential operator of the form 
D~^D, for instance, let G*-'(x) = — [a g^g^ + bSijA], where a,b > 
where A denotes the Laplace operator. The symmetric bracket (3.7) 
assume now the form 



-< ^ (f) = J (i"x|a 



V 
V- 



6^ 
(5$ 



6^ 
6^ 



(3.8) 



As we shall see this is exactly a semimetric bracket need in viscous 
fluid dynamics. □ 



We shall call semimetric dynamics a dynamics which is governed by some 
semimetric (dissipative) vector field. In local coordinate system (z^), we 
have the following system of first-order differential equations: 

=^ z\ S y= X^{z') = E G'\^)^ , i,j = 1,2,..., N, (3.9) 

i=i 

where ^ •, • ;^ is semimetric bracket, Xg is a dissipative vector field gener- 
ated by function S defined on the phase space, and G is a semimetric tensor. 
In some very special cases the function S has a physical interpretation as 
entropy. S is always non-decreasing, since 

S ='<S,S y>0. (3.10) 

If /, g are Constants for semimetric dynamics, then the Leibniz rule ensures 
that fg is also a constant, but ~< f,g y is usually not a constant due to 
the lack of the Jacobi identity. We shall now introduce a new concept a 
symmetric analogy of the Jacobi Identity which we call the SJ-identity. 
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iv) SJ-identity: 



y f,g,heT: 2^^ f,gy,hy = [^^ f,h y, g y + g,h yj y] 
^yf,heT: ^^fjy,hy = ^-<f,hy,fy. (3.11) 

A symmetric dynamical system is called a SP-dynamics ("symmetric Poisson 



dynamics") iff the symmetric Jacobi identity (3.11) holds. 



Proposition 3.2 If f,g are Constants of the SP-dynamics, then -< f,g >- 
also is a Constant of motion. In particular, if / is a constant, then ~< f, f >- 
also. 

Proof. Indeed, since -< /, 5 g,S y= we have 

^ -< f.9^=<'< f,9 ^= ^ H-< f,Sy,gy + ^-<g,Sy,fy]= Qf3.12) 

Q.E.D. 

Therefore constants of SP-dynamics form a subalgebra. This property is 
quite useful in finding constants of motion for SP-dynamics. 
Note that each autonomous dynamical system described by the system of 
the first order differential equations 

Xk = Fk{xi,. . . ,Xn) , k = l,...,n, (3.13) 

is a pseudo-metric system. Indeed, one can choose as the function S, 
S = J2k=i ^k, and as a diagonal pseudo- metric G^^ (x) = 5^^ Fi{x). 
Furthermore, locally and almost everywhere each dynamical system ( [3.13]) is 
metric. Indeed, for each point from the set {x : \/k = 1, . . . ,n : -Ffc(x) ^ 0} 
there exists a neighborhood U such that functions do not change their 
sign inside U. Let denote Sk = signF^ = ±1 in U. In the neighborhood 
U, the system should be regarded as metric system with, for instance, 
G'^{x) = Si5'^Fi{x) > and 5 = Efc 

SfcXfc. In particular, locally and al- 
most everywhere, Poisson dynamics also admit a metric description. Each 
dynamical system of the type xt = Fk{x.) should be regarded as Poisso- 
nian system after doubling the number of variables. Indeed, consider a 
canonical Poisson structure: {xi,Xj} = = {pi,Pj}, {xi,Pj} = 5ij and let 
'H(x, p) = J2kPkFki^), the canonical equations follow 

^ dF(:x.) 

i:k = {xk,n} = Fk{^), = {pj„n} = - Y,Pj-^ ■ (3.14) 
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A system is non-Poissonian (or non-metric) if it can not be written in the 
Poisson (resp. metric) form without changing the number of variables. Each 
dynamical system of the type = -ffc(x) where functions -Fa:(x) are positive- 
definite, is a metric system but it is (in the general) non-Poissonian. The 
answer to the question which Poissonian system admits a global metric de- 
scription, remains unknown. 

Morrison Q has pointed out that metric dynamical systems admit an asymp- 
totic stability at isolated maxima of function S. To show that let x be an iso- 
lated maximum of S, then certainly diS = at x, hence x is an equilibrium 
point of the semimetric dynamical system: i* =-< z^,S y = G^^ {z)djS. 
Define the function L{z) = 5(z) — 5(x), then obviously L(x) = and 
L(z) < in some neighborhood of x, since x is the isolated maximum. Fur- 
thermore, L(z) =^ L^S >- (z) =^ 5, 5 ;^ (z) > for z 7^ x and -L(x) = 0, 
therefore L is Liapunov function for the system and x is its stable equilib- 
rium point. Now if the system is metric, then L(z) S,S >- (z) > 
for z 7^ X since function S is not locally constant. Hence x is asymptoti- 
cally stable point. Generalization of the above construction to the infinite 
dimensional case is not known. 



4 Constrained metric dynamics 

In the frame work of symplectic geometry, constrained Hamiltonian dynam- 
ics can be represented by a triplet (M, A^, oj) where (M, uj) is a symplectic 
manifold, namely Phase space, and is a constraint submanifold of M. 
Antisymmetric Dirac bracket for second-class constraints |l6| is nothing 
but the Poisson bracket on some symplectic manifold A^' C A^, called the 
second-class constraint manifold |jl^ et. al. (also in [^). 
Similarly, constrained metric dynamics should be represented by a triplet 
(Af, A^, g) where g is metric tensor which is responsible for a dissipation 
and A^ is a constraint submanifold of M. We show that symmetric Dirac 
bracket for a triple (M, A^, g) is nothing but the semimetric bracket on the 
submanifold A^. It is worthy to note that any submanifold of a Riemannian 
manifold is second-class with respect to the metric bracket defined by the 
metric tensor. 

Suppose that we have a pair (M, ^) where M is a smooth manifold and ^ is 
non-degenerate symmetric (0, 2)-tensor on M. Then at each point x € M the 
map : T^M x T^M — > i? is bilinear, symmetric and non-degenerate and it 
induces a linear bijection : T*M T^M. Further, if ^ is Riemann metric 
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tensor, then is a scalar product on T^M. The non-degeneracy of the ten- 
sor ^ guarantees an existence of a (2, 0)-tensor field A : T*M x T*M — >• R. 
Let be a submanifold of M and ^|jv is supposed to be non-degenerate. 
Then at each point x of the submanifold A^, the linear space TxM decom- 
poses into direct sum of a tangent space to A^, T^N , and its orthogonal with 
respect to bilinear symmetric functional ^a;, i.e. T^M = TxN(B{TxN)^. The 
symmetric Dirac bracket with respect to a triple (M, A^, ^) is defined by 

Az)(a,/5) = CiPS,H»),PCHP)) , where P is a projection 

onto TxN along (TxN)-^, a,/3 are 1-forms on M . (4.1) 

The symmetric Dirac bracket in the space of functions is then 

^f,gyD = AD{df,dg), yf,geC'^iM). (4.2) 

If A^ is the second-class submanifold 



N = {xeM : Qi{x) = 0}, (4.3) 

then we derive explicite formula for symmetric Dirac bracket with respect 
to the triple (M, N,^). Let us denote Xg a vector field generated by the 
function g, i.e. Xg{f) = A{df,dg) =^ f,g ^= C{Xf,Xg). Let W = [Wij] = 
[-< Qi^Qj >-] = [C{Xe,,Xej)] and C = [Cij] = W'^. It is easy to see that 
the vector fields Xq. span TN^, then orthogonal projection Q onto TN-^ 
along TN has a form: 

QX = J2^iX,Xe,)CijXe^. (4.4) 

Then the orthogonal projection P onto TN along TN-^ is of the form PX = 
X — QX, hence we have PXj = Xf — J2i,j XQi)CijXQ.. Then the 
symmetric Dirac formula is of the form 

^f,gyD = AD{df,dg) = aPXf,PXg) = aXf,PXg) 

= aXf,Xg)-aXf,QXg) 

= ^ /,9 ^ ^ /,0i ^ C'ij ^ 0i,5 ^, (4.5) 



which coincides with the antisymmetric Dirac bracket formula (2.4) for Pois- 
son bracket with the antisymmetric brackets {•, •} replaced by -< •, • Our 
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procedure shown above applies to both cases, for symmetric or antisymmet- 
ric (0, 2)-tensors. 

The Dirac formula ( |4.5| ) plays a key role in the practical use of the Dirac 
brackets. Algebraically, one may use it as a definition of the Dirac bracket for 
an arbitrary symmetric or antisymmetric algebra. The disadvantage of the 
algebraic approach is based on the fact that it is very difficult to understand 
why the Jacobi identity for the new Dirac bracket holds when the above 
procedure is applied to a Poisson bracket. In the metric context, if the alge- 
braic formula (|4.5| ) is regarded as a definition of the Dirac bracket, then it is 
easy to check that new algebraic Dirac bracket is symmetric Leibniz bracket 
(i.e. the algebraic properties i), ii) hold), but the crucial non-negativity 
property iii) is not easy to verify. A simple proof will be given latter after 
the Theorem It is easy to see that 



-< ©a, / >~D= , for arbitrary function f{z) 



(4.6) 



therefore all constraints Qa ai^e Casimirs, i.e. belong to the Centrum of the 
semimetric algebra (C°°(X), -< •, • >~d)- 

Now we would like to present the new formula for calculating symmetric or 
antisymmetric Dirac bracket. 



Theorem 4.1 

The following identity holds: 



detW 



(4.7) 



where 



w 



^62,61 y 



-< Gi,GAr >- 

-< e2,@N y 



^ eAr,Gi ^ ... ^ eAr,eAr ^ _ 



The same formula holds for antisymmetric Dirac brackets. 



(4. 
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The proof is straightforward; one apphes twice the Laplace recursive formula 
for the determinant expansion to the last column and row of the matrix Wf^g . 
□ 



One consequence of ( |4.7| ) and ( |4.8| ) for semimetric bracket is that V/ we 
have -< f,f >~D= '^detvy^ ' ^"^^ therefore the inequality ^ f, f >~d^ holds, 



according to the Proposition 3.1. In one constraint case, this is equivalent 
to the Schwartz inequality. 

The theorem [4.1| provides a new effective formula for calculating Dirac 
brackets for both symmetric and antisymmetric case. Usually, the direct 
attempt to use the formula ( [4 .51 ) is impractical for a system with rather big 
number of constraints of second type. This is because it requires quite com- 
plicated evaluation of the elements of the inverse matrix Cij. Note, however, 
that when we are not interested in the Dirac brackets but only in result- 
ing equations for constrained dynamics then the evaluation of C = W^^ 
is unnecessary. Applying theorem |4.1j we find immediately equation of the 



quantity / in constrained symmetric/antisymmetric dynamics: 

/-^/,7^^z.-^-^, or f-{f^^)D-^^. (4.9) 

This formula is particularly convenient for finding constants of motion for 
constrained - symmetric and antisymmetric - dynamics. Indeed, / is a 
constant of motion for constrained dynamics iff det Wj^n = 0. 
Below we show few examples of the constrained symmetric Dirac brackets 
which are applicable in differential geometry and physics. 

Example 4.1 Consider the standard Euclidean metrics -< , >-= 5^^ = 
G*-' and the fixed surface /(z) =0 in R". Here the function / is assumed 
to be smooth and with zero as it regular value, that is /~^(0) is a close 
regular n — 1 dimensional differential submanifold in R" . Using /(z) = as 
a constraint we find the Dirac semimetric brackets 

^z\z^ '^D=5'^ -n'n^ , (4.10) 

where n(z) = is an unit normal vector to the surface at z. The metric 

tensor G^(z) = 5"^^ — n^n^ is nothing but the induced metric tensor on this 
surface. □ 

The following example illustrates how to derive the metric structure for 
lattice-spins that is a set of classical spins Sa where a labels the lattice 
sites 
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Example 4.2 We introduce the usual lattice-spin metric brackets as 

^SlSiy = 6abS'^\Sa\ = G'^„ i,j = 1,2,3, a,b = 1,2,..., NiA.ll) 
and we define 2A^-dim surface V by the following system of N constraints 



G„(S) = \Sa\' -rl = J2iSlf -rl = 0, a=l,2, 

i=l 

The Dirac metric brackets for the surface V are 



.N. 



SI 



(4.12) 



(4.13) 



□ 



The next example illustrates the metric structure for energy-conserving rigid 
body. 

Example 4.3 

One may postulate the metric brackets for a rigid body as 



5ijK{uj) , where K is some function of the rigid body 



angular frequency uj . 
a) Consider energy as constrained surface 

3 

k=l 

Calculating the Dirac metric brackets one gets 



-< uji,ujj y{e}= K{uj) 



We can consider some particular cases: 



liljtditdj 



(4.14) 



(4.15) 



(4.16) 



al) One may choose K = X)l=i -^I'^fc) where Ik are moments of 
inertia with respect to main axes of the rigid body. Then the 
Dirac metric brackets are 



-< tdi,tdj ;^{0}— <5f 



k^k 



liljtditdj 



(4.17) 



This metric structure coincides with the metric structure postu- 
lated by Morrison 0]. 
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a2) Since the Poisson structure of rigid body is the same as for clas 
sical spins, we may postulate K = \uj\, where \uj\ - 



The Dirac metric brackets fohow 



-< ^{e}= I'^l 



liljLUiUlj 



J2 k=l H 



OJt 



(4.18) 



b) Consider Poissonian Casimir as constrained surface 

k=l 



The Dirac metric brackets fohow 

^ uji,ujj ^{e}= K{uj) 
For instance, K{u)) = \io\, we get 



IjOiUJj 



OJiUJ 



(4.19) 



(4.20) 



(4.21) 



□ 



Next examples show how our formalizm works in the Hilbert spaces. 

Example 4.4 Let Ph = L2(R";R^) be a Hilbert space of real, vector val- 
ued square integrable functions with standard scalar product < • | • >, i.e. 
< f I g >= /d"x f(x) • g(x). where f = (/i, ...,fd), g = {gi, ■ ■ ■,gd)- Let 
||f|P =< f I f >. In the space of all smooth functionals over Ph, according 



to the construction given in the Example 3.2, the semimetric structure may 
be defined by 



-< 91,92 >~ 



(f ) = x: / rf"x 

i=l 



<5/i(x) <5/i(x) 



h,(l)2-Ph^R, (4.22) 



where denotes a Gateaux functional derivative. 

o/fe 

Consider a surface of infinite dimensional sphere S°° with radius 
subspace with one constraint 



5°° = {f e Ph 



(4.23) 
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Now we calculate the Dirac metric structure for the sphere S°°. The metric 
bracket ( [4.22 ) can be rewritten, introduce the canonical metric tensor G 



G^^(x,y) =^ /.(x),/,(y) ^= 5,,<5(x - y) , 



(4.24) 



as 



X(/<i,</.2 ^ (f) = Jd^xd^y^ 



SM^)6fjiy) <5/,(y)<5/,(x) 



The Dirac semimetric structure on follows 



^ /i(x),/j(y) = %5(x-y) 



/»(x)/,(y) 
llflP 



-< /*(x),/i(y) >- . 

(4.25) 



(4.26) 



□ 



Example 4.5 Let ^' = ^'i + i^2 and its complex conjugate = — i"^2 
where ^i, ^2 are real functions integrated by square, i.e. they belong to the 
Hilbert space L^. We define metric structure by 

^^k{^),^i{y)y=^6kiSix-y), where fc,/ = 1,2. (4.27) 
One can rewrite it to the form 

^ ^(x), M/(y) vl/*(x), M/*(y) y= 0, ^ M'(x), M/*(y) y= 5(x - yX4.28) 

which we call by canonical metric bracket for Quantum Mechanics. The 
Dirac structure on the sphere ||^|| = Const in the Hilbert space follows 



2||^'| 



^ M/(x), M/*(y) ^Z5= <5(x - y) - ^-^^^4.29) 
The Dirac brackets of the components one get in the accordance with 



^ ^'fc(x),^';(y) ^D-- 



6ki6{x-y) 



^fc(x)^;(y) 

ll^'IP 



(4.30) 
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□ 



Example 4.6 (The canonical description for incompressible, viscous fluid 
dynamics is based on this example) 

Let Ph = W^(i'2)(R"; R") be a Sobolev space of real functions. In the space 
of all smooth functionals over Ph we introduce a semimetric structure 



-< 4>i,(l>2 (J) = / d"x<^ a 



(5J(x) 



V 



<5J(x) 



^1 



V 



M;(x) 



where J denotes the real vector in the n-dimensional Euclidean space, and 
a, b are real non-negative coefficients. We can rewrite this semimetric struc- 
ture in the form 



-< J,(x), Jj(y) y 



OXi OXj 



5(x - y) . (4.32) 



Consider infinite-dimensional subspace V of divergence-free functions (in- 
compressibility condition) as a system with infinite number of constraints: 



V = {J ePh : Vx ex(J) = Vx • J(x) = 0} . 
Then the Dirac semimetric structure for the subspace V follows: 

d d 



(4.33) 



-< Ji(x), Jj{y) )~D 



dxi dxj 



5(x-y) 



J dzdz' ^ j,(x),e(z) y c{z,z') e(z'), Jj(y) y, 

(4.34) 

where C is an inverse symmetric operator of the constraint matrix 

C(x,y) = ^ /dzG(|x-z|)G(Iz-y|), (4.35) 
a + J 

here G(|x — y|) denotes the standard Green function (fundamental distri- 
bution) of the Laplace equation, i.e. AxG(|x — y|) = S{x — y). Putting 
eq. ( [4. 351 ) back to eq. ( |4.34| ) and after some simple calculations finally we 
obtain 



-< Ji(x), Jj{y) yn 



52 

(Jjj Ax - ^ 7] 

OXiOXj 



5(x - y) . (4.36) 



(4.31) 
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Physically, the last eq. ( 4.36 ) fully describes dissipative structure for incom- 
pressible viscous fluid. We shall see that more clearly in subsection |6.5| . 



□ 

5 Semimetric-Poissonian systems 

The physical systems are usually dissipative. It turns out that both Pois- 
sonian and metric structures alone are not enough to describe dissipative 
systems. However, a proper combination of these two types of dynamics 
can, for many interesting cases, provide a satisfactory and fully algebraic 
description of dissipative dynamics. 

A semimetric-Poissonian bracket on the space of functions J- = C°°{X) 
is a bilinear operation {{■^■}} : T x T ^ T which is a linear combination of 
a Poisson and a semimetric bracket 

V/, <? e : {{/,<?}} = ^ /,5 ^ , (5.1) 

where {•, •} is the Poisson bracket and ^ ■, • is the semimetric bracket. 

Definition 5.1 Semimetric-Poisson manifold is a pair (Af , 11 — G) where 11 
is a Poisson tensor, G is a semimetric tensor. 

We shall call a Semimetric-Poissonian dynamics a dynamics governed by 
the following system of equations: 

= {{z\ CD}} = {z\ $}- -< z\ ^ y= X^{z') - Xiiz^) , (5.2) 

where ^> is some phase space function. In real physical applications it is 
often the case that $ has the interpretation of the system free energy. It is 
then a matter of convention to chose minus sign in the eq.( |5.2D . Indeed we 
have then 

6 = {$,$}- ^ ^= ^< 0. (5.3) 
what describes the dissipation of energy. 

It is often even more convenient to go a step further and decompose function 
^> into two parts: the internal energy Ti and the dissipation function S, so 
^ = Ti. — S. Hence, when 5 is a Casimir of the Poissonian part, {/, 5} = 
for all /, the evolution of some "observable" / in the semimetric-Poissonian 
dynamics follows 

/ = {{f,'^>}} = {f,n}-^f,n-sy. (5.4) 
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Equation ( |5.4D for f = S give us the semimetric-Poissonian formulation for 
the "second law of thermodynamics" , namely 

S+^S,ny = ^S,Sy>0, (5.5) 

where the expression on the l.h.s. of Eq.( |5.5| ) is just the convective time 
derivative of S along the time trajectory in the semimetric-Poissonian phase 
space. 

Note that if x is an isolated minimum of <I>, then the function L{z) = 
^{z) — $(x) is a Liapunov function for semimetric-Poisson system. Hence, 
we obtain 

Proposition 5.1 If x is an isolated minimum of the free energy function 
then x is a stable equilibrium point for semimetric-Poissonian system 
z = {z,<^>}— ~< z, $ y. Furthermore, if the system is metric-Poissonian, 
then X is asymptotically stable point. 



Example 5.1 Consider a modification of the harmonic oscillator described 
by 

xi = X2 — axi{x1 + X2), X2 = —xi — ax2{xi + . (5.6) 
This system is semimetric-Poissonian with 



2 I 2 
X ^ ~t~ X 2 



{XI,X2} = 1, $ = ^ 

-< xi,xi y= axf, -< X2,X2 y= axl, ~< xi, X2 y= axiX2 ■ (5.7) 



Point (0, 0) is an isolated minimum of <5, hence it is a stable equilibrium 
point for the system. □ 

The proposed in preceding sections formal structure of Poissonian and semi- 
metric dynamics allows us to suggest the following scheme for construction 
of constrained dissipative dynamics of real physical systems. 

a) Consider a canonical Poisson structure and semimetric structure. The 

semimetric structure must be postulated according to our physical 
insight in the nature of relevant dissipative processes. 

b) Choose sets A, B of second-class constraints for the Poissonian and 
semimetric part, respectively. Note that any constraint is second-class 
for non-degenerate semimetric structure. 
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c) Calculate antisymmetric and symmetric Dirac brackets with respect to 
the set of constraints A and B and later take a union of them. 

Most interesting cases take place when the set of constraints for the semi- 
metric part is a subset of constants of motion for the Poissonian part, i.e. 
system is dissipative however some of constants of the Poisson dynamics re- 
main constants for the semimetric-Poissonian dynamics. One may use this 
feature to design many interesting dissipative systems. As an example, we 
illustrate how to design variety of damped rigid body dynamics in the next 
section. 

6 Applications: Physical Examples 

6.1 Pcirticle on the hypersurface with friction 

We shall begin our illustration of the Dirac bracket applications by discussion 
of a simple example, namely the classical particle moving with friction on a 
hypersurface iSconfig = {x € R"^ | /(x) = 0}. Denote the particle position 
by X and its conjugate momentum by p and let the friction force experienced 
by that particle be proportional to the particle velocity. The phase space 
R^" is now equipped with two structures: canonical Poisson and semimetric 
structure. Then the Poisson structure follows is described by 

{xi,Xj} = = {pi,Pj}, {xi,Pj} = Sij , (6.1) 

and the semimetric structure is defined by 

~< Xi,Xj y= =-< Xi,pj y, -< pi,pj y= 6ijXi , (6.2) 

where Ai(x) > is the directional and space dependent damping coefficient. 
The semimetric- Poisson structure is defined by {{•,•}} = {•,•} < •,• >-■ 

2 

Consider W = ^ + ^(x), the dissipative dynamics derived from the above 
structures follows 

Pi Pi 

Xi = {{xi,n}} = —, Pi = {{Pi,'H}} = -Ai(x)— - — , (6.3) 
m m oxi 

which can be rewritten in the Newtonian form mxj -|- Aj(x)i;j — Fj(x) = 0, 
where F(x) = — ^ is a potential force. In the particular, when Ai(x) = 
A(x), we have mx + A(x)x — F(x) = 0. 

Next consider a fixed algebraic surface /(x) = in R". We make assump- 
tions that / is smooth and zero is its regular value. The second assumption 
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ensures that <5gQQgg = /~^(0) is a close regular n — 1 dimensional differential 
submanifold in R". Moreover this assumption guarantees V/ / 0, so we 
can use the gradient to define the normal vector on >5(,Qj^gg. 
The set of constraints consists now of two elements: 



Gi = /(x)=0, 02 = p- 1^ = 0. 

ax 



(6.4) 



For the Poissonian dynamics both the above constraints Qi are second-class 
in the Dirac classification. 

Denoting the unit normal vector to the surface / at the point x by n(x), 
n(x) = TSTjl^ , the antisymmetric Dirac brackets for the Poissonian part 

I I 

of our construction are 

{xi,Xj}D = 0, {xi,pj}D = Sij ^^-^ = Sij -ni{x)nj{jc), 



91 
dx 



{Pi,Pj}D = 



r d 1 


dl 


_^ 




d ■ 


9/1 


p ^ 

ox. 


dxi 


dxi 


P 


9x. 


dxjj 



p • 



d_ 



ni X 



nj(x) 



P • 



d_ 
5x 



(6.5) 



For the semimetric dynamics only the second constraint B2 is second-class. 
The symmetric Dirac brackets for the semimetric part 



-< Xi, Xj >D= =< Xi,Pj , -< Pi,Pj ^D= SijXi - 



Aj(x)Aj(x)^^ 



Efc Afe(x) 



dxk 



(6.6) 



In particular, when Ai(x) = A(x), the symmetric Dirac brackets can be 
written in the form 

-< Xi,Xj >-D= =< Xi,pj >-D , -< Pi,Pj ^D= A(x)[(5ij - ni(x)nj(x)] . (6.7) 

Finally take an union of these two Dirac structures 

{{xi,Xj}}D = 0, {{xi,pj}}D = 5ij - ni{x)nj{x) , 

{{Pi,Pj}}D = 



1 9/ 




d ■ 


dl 






d ■ 


5/1 


\dxj 


P 


5x. 


dxi 


dxi 


P 


5x. 


dxjj 



-A^(x) 



Eifc Afe(x) 



ML 

dxk 



(6.^ 
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When the Hamiltonian for that system has the form = ^ + ^(x), the 
dissipative Hamilton-Dirac equations of motion follow as 



Xi 



Pi 



1 Pi 

{xi,n}- -< Xi,n >-= — \pi- {p- n)ni] = — 
m m 



{Pi,n}- ^ Pi,n y=Fi 



1 

F • n H p 



m 



P 



d_ 
9x 



n 



Ai(x) 



m 



df I Y.jH^)U-Pj 



Pi 



dxi 



dxi, 



For isotropic damping, when Ai(x) = A(x), we can rewrite these equations 
of motion in the Newtonian form as 



m X + A(x)x 



F(x) 



F(x) • n(x) + mx • — n(x) 



n(x). (6.10) 



6.2 Variety dynamics of damped rigid body 

The usual Poisson brackets for a rigid body angular velocity vector cD, after 
suitable rescaling, are 



k 



(6.11) 



Suppose that there is no second-class constraints for the Poisson part of 
the rigid body dynamics and just one such a constraint for its symmetric 
part. Assume now that this constraint is such that the system energy (or 
any Casimir function for the metric part) is constant. The choice of that 
constraint determines details of the semimetric-Poissonian structure. 
Using the canonical metric brackets (4.16) and the constrained energy Ti 
from example 4.3 we obtain Dirac metric bracket for rigid body (4.17). 
Combining these two structures one finds the metric- Poissonian brackets for 
rigid body M: 



{{uJi,ujj}} = {uJi,ujj} - X ^ LLii,ujj yo 

/ 3 

= ^ijk^k — A 5; 

I \k=l 

Consider the system "free energy" 



liljUJiUJj 



(6.12) 



^=n-s=\ 



.k=l 



■5(kP), 



(6.13) 
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The equations of, energy conserving, motion for the damped rigid body are 

= u;2W3(/2 - -^3) + 2A5'a;i hih - h)'-^l + h{h - h)^l 

and its cychc permutation . (6.14) 
The dissipation of the free energy follows 

6 = -S = -X<S,S '^D 
dS 2 



-A 



duj 



\k=l / 



liOJ, 



95 \ 2- 

dui) 



<0, 



(6.15) 



and there is no dissipation of energy iff Ji = /2 = -^3 • ^ 
Alternatively, we can derive new semimetric-Poissonian bracket by combin- 
ing metric bracket ( [4.18| ) with the standard Poisson bracket for the rigid 
body 



{{uJi,UJj}} 



{uJi.UJj} - A -< UJi,LOj >- 



D 



^ijk^^k — A|a;| 



liljUJiUJj 

t2, ,2 
Z^fc=l ^k^k 



(6.16) 



Equations of motion following from the free energy ( |6.13| ) are 

dJi = {{uji,<^}} = {uji,n} + x^uji,s 

h{l2-h)0jl+h{,h-Il)ul 



^2^3 (-^2 - -^3) + 2\S'\uj\uJi- 



l^k=l ^k^k 



and its cyclic permutation . □ (6-17) 

Finally we derive new semimetric-Poissonian bracket by combining metric 
bracket (4.21) with the standard Poisson bracket for the rigid body 

{{uJi,UJj}} = {uJi,UJj} - X UJi,UJj 



^ijk^k — A|u;| 



Si 



UJiljOj 



a; 



(6.18) 



Again using the free energy ( |6.13| ) we find 

cji = {{uji,^}} = {uji,n}-\<uji,n>-D 



X\uj\loi 



h 



and its cyclic permutation 



(6.19) 
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The above dynamics does not conserve the energy, but |u;p remains the 
system Casimir. The energy dissipation is given by 



-A ^ n,nyD-- 



-X\uj\ 



2 ,2 {T^kh^ff 



k^k 



\UJ\ 



<0, (6.20) 



since i^k^k^Vl I'^P — {Ylk-^k^'kf ■ For A > 0, \uj\ > 0, there is no dissipa- 
tion of energy when Ii = I2 = I3. □ 



6.3 Classical damped spins 

The Poisson bracket for lattice-spins 



^ab^ijkS(i , 



(6.21) 



where a labels the spin location and k = 1,2,3, guarantees that length of 
each spin |Sa[ is a Casimir. Our model is based on the assumption that 
there are no second-class constraints for the Poissonian part and only one 
constraint for the metric part. 

Starting from canonical metric bracket ( 4.11| ) and using the length of the 
spins I Sal constraint, as in the example |4.2| , we obtain the Dirac metric 
bracket for lattice-spins ( |4.13 ). Combining these two structures one finds 
the metric- Poissonian of lattice spins [S] : 



{{SI si}} = {SlSi}-X^SlSiyD 



^ab^ijkSa — AJ*-' |Sa| 



(6.22) 



'->a'~>b 
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where A is the damping coefficient. 

Now we can easily derive the Landau-Lifshitz-Gilbert equation of classical 
damped lattice-spins 



Sa 



{{Sa,n}} = SaXB,f, 



A 



Sg X (Sg X Be/,a) 
|Sa| 



(6.23) 



where Bef,a = ~'§^ is the effective magnetic field acting on the spin Sg. 
The dissipation of energy follows 



n 



-A -< Ti.,Ti. >~D= ~A ISgl 



■d2 (-^e/,a • Sg) 
^ef,a 



|Sa| 



<0. (6.24) 
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6.4 Dissipative quantum mechanics 

Combining canonical Poisson brackets for quantum mechanics [QM] 

{vl/(x), ^(y)} = {M'*(x), **(y)} = 0, {^(x), vl.*(y)} = lj(x - y) , (6.25) 

and constrained metric brackets for Q M ([4.2g| ), which were derived from 
the canonical metric bracket for QM ( 4.28| ) for the physically important 
constraint, namely the conservation of the wave function norm (probability 
conservation) as in the example 4^, we find dissipative Metric-Poissonian 
structure for QM : 



{{^(x),*(y)}} 
{{^*(x),^*(y)}} 
{{vl/(x),^*(y)}} 



A A'I'(x)^'(v) 
{^(x), ^(y)} - - ^ vl/(x), ^{y) yD= ^ 2\m' ' 

{^I/*(x),^*(y)} - ^ ^ ^*(x),^*(y) yD= ^^^p^ 



— <5 x-y - - 

in n 



5(x-y) 



^(x)^*(y) 



2||^'| 



(6.26) 



where A is the (undefined) damping constant. 

Using the conventional Hamiltonian for the Schrodinger equation 

7^(*,^*) =< >= J d"x^*(x)i?\^(x), (6.27) 

where H is the quantum mechanical (self-adjoined) operator we obtain the 
evolution of the wave function in the form: 



ih{{^{x),n}} = H^{x) + iX 



II^IP 



H 



^'(^)28) 



This equation is known as the Gisin dissipative wave equation |11]. Here the 
construction of the semimetric bracket ensures that the norm of the state 
vector is reserved (so the probability is conserved) since it is a second-class 
constraint for the semimetric part. The dissipation of energy follows 



n = - ^n,ny=^ 

n 



iH'^lr + fTTTT^ 

hi' 



< 0. 



(Schwartz inequality) 



(6.29) 
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In the above equation the equality is achieved for ^' which are the eigenstates 
of the Hamihonian. The damping in the Gisin equation refers, therefore, to 
the transition ampUtudes only. When the initial wave packet is constructed 
from eigenstates corresponding to energies E > Eq then the final state of 
the evolution described by eq. (6.28) is the eigenstate with energy £"0. 
This property distinguishes the Gisin dissipative wave equation from other 
dissipative Schrodinger equations. 



6.5 Viscous fluid dynamics 



In the fluid mechanics the state of an isothermal fluid is described by its 
mass density and velocity fields (^J, u) or by (i?, J) where the current J field 
equals J = qw.. The semimetric structure for fluid dynamics pi follows 



-< Jfc(x), Ji(y) >~ 



^ ^(x),e(y) ^= 0, ^ ^(x), Jfc(y) ^= 0, 
92 



dxkdxi 



(6.30) 



where a = C + 3 and V are bulk and shear viscosity respectively. 
The structure ( |6.30D is semimetric, indeed \tQ we have: 



k,l 
fix < a 



Q2 



V 



60 



dxkdxi 
V 



k 



60 



5(x-y) 



6g 



5Ja:(x) 



6Ji{y) 



(6.31) 



Note, that this is not a metric bracket, since any functional which depends 
only on g, has zero semimetric bracket with himself. Also note that the 
kinetic fluid energy E^in = J dx J'^/2g dissipation follows directly from 



-^kin 



^ -^kini -^kin ^ 



dx< a 



V 



q{x) 



V 



q{x) 



< 0. 



(6.32) 



In classical hydrodynamics a particular role is played by the incompressible 
fluid assumption, for incompressible viscous fluid according to eq. ( 4.36| ) we 
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have 



^ ^(x),^(y) 0, ^ Q{x),Ji{y) ^0=0, 



^ Ji(x), Jj(y) >-D= -rj 



Sij^x - ^ ^ 

OXiOXj 



<5(x-y) 



(6.33) 



The dissipative energy for viscous incompressible fluid should be easily cal- 
culated 



Eki 



duk ^ duj_^ 2 
dxi dxk , 



< 0. 



(6.34) 



Finally, put ^ = the kinematic viscosity, using (6.33) and the Poisson 
part calculated in |^ we easily derive non-local evolutional equations for 
incompressible viscous fluid 



'dt 



+ (u • V)u - z^Au = V j dz G(x - z)Vz [(u • V)u - z^A^u] 



(6.35) 



7 Final comments 



The description of the dissipative systems dynamics is usually undertaken 
within the framework of non-equilibrium statistical mechanics. The "simpli- 
fied" version of the full many-body description often used in applications, for 
example in phase transformations physics, statistical theory of turbulence, 
granular media dynamics etc. is the kinetic equation for the time evolution 
of the "relevant degrees of freedom" distribution function. This equation is 
derived making strict assumption about the nature of the fluctuations in the 
system - the underlying stochastic process performed by the relevant sys- 
tem degrees of freedom in the full phase space of the physical model. The 
equations of motion for the coarse grained variables describing the meso- or 
macroscopic system properties are obtained from that kinetic equation by 
one of the known procedures, which once are well justified and understood, 
and occasionally are just a heuristic chain of semi-mathematical operations. 

In this paper we have discussed a novel approach to description of the 
dissipative systems dynamics, which is purely algebraic. Instead of making 
series of assumptions on the level of microscopic physics of the problem we 
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assume existence of certain algebraic structure, akin to that used in Hamil- 
tonian formulation of classical dynamics which permit us to derive meso- or 
macroscopic dissipative equations directly from the system free energy. The 
basic ingredient of that procedure, the semimetric Leibniz bracket for dy- 
namics variables over the whole phase space, is postulated according to our 
knowledge about dissipative processes. With knowledge of the Poissonian 
structure for non-dissipative part of the system dynamics, the symmetries 
of the problem and using Dirac machinery [l[ one can derive easily different 
Dirac structures, which are otherwise hard to postulated, describing dissipa- 
tive dynamics. One may use the algorithm proposed here to design variety 
dissipative dynamical systems with required conservative observables. 

As shown in the paper this permits us to build up the metriplectic dy- 
namics scenario for several non-trivial systems: classical particle physics, 
many spins dynamics, rigid body dynamics, compressible and incompress- 
ible viscous fluid dynamics and some quantum mechanical problems. Several 
other applications, notably the relativistic charged particle systems which 
can also be formulated within the metriplectic scenario have been discussed 
previously [0, The quantum applications are of particular interest in 
view of some similarity between the metriplectic approach and the Lindblads 
construction ||l^ - the standard tool in dissipative quantum system analysis. 
We expect to comment on the connection between both these approaches in 
the forthcoming publication. 



Appendix 



A symmetric Leibniz bracket which satisfies the SJ-identity is called the SP- 
bracket. The tensor G which generates SJ-bracket is called SP-tensor. In 
the local coordinates 



V, 



ij=l,2,...,Af 



N 

E 

k=l 



Q^k Q-^k 



(7.1) 



Hence, SP-bracket and SP-tensor are symmetric analogy to Poisson bracket 
and Poisson tensor, respectively. In particular each symmetric tensor where 
Gij are constants (i.e. do not depends on points), generates a SP-bracket. 



For a non-trivial example of SJ-tensor, see Example 7.1. 

SP-manifold is a pair (M, G) where G is a symmetric SP-tensor of the type 

(2,0). 
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Example 7.1 It is easy to sec that G = Ylii,j=i ^i^j'^ ® ^ is SP-tensor. 
Hence, (R^,G) is SP-manifold. □ 

Similarly, with the concept of SJ-identity one can define SL-algebra, SP- 
algebra which are symmetric analogy of Lie algebra and Poisson algebra, 
respectively. To prove that there is no another symmetric Jacobi identity 
we need few elementary algebraic concepts. By an identity of the algebra !F 
we mean a polynomial P in some free algebra which is identically zero when 
the generators are replaced by any elements of We are interested on the 
3-linear identities in which each term involves two pairs of brackets, i.e. of 
the form: a f,g y,h y +b f,h y,g >- +c -<-< g,hy,f y= 0. Now 
we would like prove the following result: 

Proposition 7.1 There arc exactly two types of 3-linear identities for sym- 
metric algebra in which each term involves two pairs of brackets: 

a) The SE-identity (symmetric version of the Engel identity): V/ G J^: 

b) The SJ-identity (symmetric version of the Jacobi identity) : 

V/,5: ■<-^ f,f y,9y=<^ f,9yjy^ 
yf,g,h: 2 f,gy,hy= [-^^ f,h y,g y + g,h y, f . 

Proof. Suppose that a -<~< f,gy,hy+b -<~< f,hy,gy+c ~<~< g,h >- 
, f y= 0. For f = g = h we have (a + b + c) -<-< f,f y, f y= 0, therefore 
the identity must be of the type a) or a + 6 + c = 0. If a + 6 + c = 0, then 
for f = hwc have (a + c) f, g y, f y +b -<~< f,fy,g y= 0, i.e. it must 
be of the type b), since a + c = —b. 

Q.E.D. 

Similarly, one can prove that there are exactly two 3-lincar identities in 
which each term involves two pairs of brackets for antisymmetric algebra: 

a) The Engel identity: 

^f,9-- {{f,9}J} = ^yf,g,h: {{f,g},h} + {{h,g},f} = 0. 

b) The Jacobi identity: 

V/, 9, h: {{f,9},h} + {{g, h}, f} + {{h, f}, g} = Q . D 
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Even the Jacobi identity and SJ-identity are formally similar as we have 
seen above, but there exists also a fundamental difference between them. 
Indeed, the former is equivalent to Xh{f,g} = {Xhf,g} + {f,Xhg}, i.e. 
Hamiltonian vector fields act as derivations for Poisson bracket, but for the 
latter ^ f,g Xj^ f,g + ^ f^^h9 the other words, the 

Jacobi identity manifests some basic features like cyclicity, derivation and 
signature property while the S'J-identity does not. 

Furthermore, the most natural symmetric bracket ~< A,B y= 
does not satisfy the SJ-identity in the general. Hence, there is no natural 
representation for 5P-algebras. 
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